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ABSTRACT 

In this paper, we review the theories and observations which together have led to the 
concept of the Baryon Catastrophe: observations of the baryon fraction on the scale of 
clusters of galaxies appear to be at least three times as high as the universal baryon 
fraction predicted by the theory of primordial nucleosynthesis in a flat, SIq = 1, universe. 
We have investigated whether this discrepancy could be eliminated by treating the 
intracluster g multiphase medium, and found that this treatment both lowers 

the calculated mass of gas in a cluster and increases the inferred gravitational potential. 
These combined effects can reduce the calculated baryon fraction by between a quarter 
and a half: the precise amount depends upon the volume fraction distribution of density 
phases in the gas but is independent of the temperature profile across the cluster. Thus 
moving to a multiphase intracluster medium cannot resolve the Baryon Catastrophe by 
itself; other possible causes and explanations are discussed. 
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lOl INTRODUCTION 

0^ 

— I The concept of the Baryon Catastrophe was first formulated 
when it was realized that clusters of galaxies contain more 
|ihan their fair share of baryonic matter. Primordial nucle- 

Q osynthesis is one of the most well understood theories of 
modern physics and predicts that the density of baryonic 
matter in units of the critical closure density, Qb, lies in the 
range, 0.010 < Ubh^ < 0.015 (Walker et al. 1991), where 
Ho = 100/tkms~^ Mpc~^ . The theory of inflation predicts 
that the density of the universe is very close to the critical 
density, Q = 1, and this is supported by observations of pe- 
culiar velocity fields and large scale structure (see e.g. Dekel 
1994). If these two theories are to be compatible with one 
another, then non-baryonic matter is needed in addition to 
the baryonic matter, with a density of ~ 100 — 1.25h~^ per- 
cent of the critical density. Even if a very low value of the 
Hubble parameter is assumed, such as A = 0.3, the Universe 
should still contain only 14 percent baryonic matter com- 
pared with 86 percent non-baryonic matter, giving a baryon 
fraction of 0.16. 

The largest scale on which the baryon fraction can 
be reliably calculated is that of clusters of galaxies. Re- 
cent observations (Fabian 1991; White et al. 1993; Bohringer 
1994; Henriksen & Mamon 1994; Allen et al. 1995; Buote & 
Canizares 1995; Elbaz, Arnaud & Bohringer 1995; White & 
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Fabian 1995) have found that clusters of galaxies contain 
three times as many baryons as predicted by the theory of 
primordial nucleosynthesis in a critical density universe. 

We are left with the following options: (i) the observa- 
tions, or their interpretation, are incorrect, (ii) Qo < 1, or 
(iii) the limits from primordial nucleosynthesis must be re- 
laxed (see e.g. Sasselov & Goldwirth 1995). This paper tests 
the first of these possibilities. We show that treating the in- 
tracluster medium (icm) as multiphase (in which there is an 
emulsion of gas blobs of differing density at any given radius) 
can lower the estimate of the baryon fraction by as much as 
26 to 49 percent. The following section provides the moti- 
vation for and a formal description of the multiphase icm. 
Our cluster model is described in Section 3 and the results 
presented in Section 4. The implications for, and possible 
resolutions of, the Baryon Catastrophe are discussed in the 
final section. 



2 THE MULTIPHASE INTRACLUSTER 
MEDIUM 

2.1 Motivation: cooling flows in clusters of 
galaxies 

The treatment of the intracluster g multiphase 
medium is motivated by the study of cooling flows in clus- 
ters of galaxies. The central gas density in clusters is ex- 
tremely high, therefore close encounters between particles 
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are frequent, and energy is lost due to the emission of X-rays 
by thermal bremsstrahlung and, at lower temperatures, by 
line radiation. The cooling process is thermally unstable and 
overdense blobs of gas will rapidly cool further, condensing 
into molecular clouds or low-mass stars. As matter in the 
centre of the cluster cools, more matter moves inwards to 
maintain pressure support, so forming a cooling flow. For a 
review of cooling flows see Fabian (f994). 

The amount of gas deposited in these cooling flows 
is considerable, up to f 000 Mq yr~^ , (Thomas, Fabian & 
Nulsen f 987), and may contribute to the formation of a cen- 
tral dominant galaxy. The matter is deposited not solely in 
the centre of the cluster, however, but over a wide range of 
radii, out to at least f50/t~^kpc. In order for this to hap- 
pen, there must be a range of densities and temperatures at 
each radius in the cluster. If this were not the case and, as is 
usually assumed in cluster models, there were a single den- 
sity and temperature corresponding to each radius, then the 
gas would cool catastrophically at the centre of the cluster, 
giving a much more centrally-condensed surface brightness 
profile in X-rays than is observed. 

The cooling flow models only constrain the structure 
of the icm within the region in which the cooling time is 
much less than the age of the cluster. There is no reason to 
suppose, however, that a single-phase model is appropriate 
at larger radii. The history of the icm is poorly understood 
but is likely to be very complex. In particular, the moderate 
metallicity (approximately half solar) suggests that it has 
been lost from galaxies either by stripping or by galactic 
winds. In addition, mergers of small groups and subclusters 
may introduce and mix together gas of a wide range of en- 
tropies. If the gas phases were free to move past one another 
then convection would stratify them on a dynamical time. 
However, the cooling flow observations show that this does 
not happen, at least in the cluster core. There are a variety 
of physical processes which may prevent convection such as 
magnetic pinning, or simply viscous drag (see the discussion 
in Nulsen f986). Here we assume that the icm is multiphase 
throughout in order to see what effect this has on estimates 
of the baryon fraction. 



2.2 Description: the volume fraction 

The formalism for describing a multiphase icm was set out 
by Nulsen (f986). He introduced the volume fraction distri- 
bution, f{p,r,t), such that fdf denotes the fraction of the 
volume at radius r and time t which contains gas in the den- 
sity range [p, p + dp]. Under the assumption that the density 
phases at any particular radius are in pressure equilibrium 
(which will be a good approximation until the cooling time 
drops below the sound-crossing time of the individual gas 
blobs) and that they are comoving, then it is possible to 
derive multiphase versions of the fluid equations which de- 
scribe both the motion of gas through the cluster and the 
time-evolution of the volume fraction. Throughout the bulk 
of the cluster, where the mean cooling time is much longer 
than the age of the cluster, the distribution of / will change 
very slowly and, in the absence of a detailed model for the 
formation of the icm, is poorly constrained. Within the cool- 
ing flow, however, / is strongly modulated by cooling such 
that, at high densities. 



(1) 



where a is the slope of the cooling function (see below). 

A heuristic argument for the above form of the distribu- 
tion function is as follows. Consider an initial distribution, 
/j. If the phases of material at the highest densities of 
are considered, then their cooling time is relatively short. If 
this timescale is less than a Hubble time, then these phases 
will already have cooled, and will not feature in the volume 
fraction distribution function today. On the other hand, the 
phases of material at the lowest densities will have a cool- 
ing time which exceeds the Hubble time, and will not have 
cooled by the present. We are therefore left with a very nar- 
row range of initial densities which are in the process of 
cooling today, for which the value of is approximately 
constant. The functional form of the highest density phases 
will therefore be approximately independent of the initial 
distribution. 

The evolution of / is governed by the Energy Equation 



dt p 



(2) 



where T is the temperature, s the specific entropy, and { 
is the emissivity (i.e. the emission per unit volume). For an 
Equation of State 



P : 



(3) 



where P is the pressure and u the specific thermal energy of 
the gas, then 



3 A; , / P 
s = In 



2 fimj 



5/3 



(4) 



where k is the Boltzmann constant and finiH the mass per 
particle. Writing 



{ = ra^A(T) = ra^AoT" 



(5) 



where n = p/fimH and « = §■ for thermal bremsstrahlung, 
then a little rearrangement gives 



P _3P 2 ra^A 
p~ 5 P~^ 5 P ' 



(6) 



At high densities the second term on the right-hand-side of 
the above equation dominates. Then 

P-g— . (7) 

where we have used the fact that the phases are in pressure 
equilibrium, pT =constant. 

The equation of mass conservation in spatial and den- 
sity space is 

§-^ipf) + ^-ipfv) + j-^ipfp) = 0. (8) 

Unless the functional form of / changes very rapidly (which 
we have argued above is very unlikely) then the final term 
dominates at high density, giving 



pfp ~ constant. 



(9) 



or / ~ as in Equation f . 

The above argument is only suggestive but was given 
some credence by the simulations of Thomas (f988) who 
followed numerically the time-evolution of a range of volume 
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fractions. In every case a high-density tail of the correct form 
was reproduced. In fact, all volume fractions, /j, which were 
initially restricted to a narrow range of densities, rapidly 
converged towards a pure power law. 



fi 



(3- 




-(4-a) 



P > Prah 
P < Prah 



(10) 



where the constant of proportionality has been chosen so as 
to give a total volume fraction of unity when integrated over 
all densities. A multiphase, constant-pressure cooling flow 
model with this volume fraction distribution gives a mass- 
deposition profile M (X r^^* which is conveniently close to 
that expected for the distribution of dark matter around the 
central dominant galaxies in clusters. This suggests that fi 
may be ubiquitous in nature and may be taken as a canonical 
form for the volume fraction. However fi is just one of a 
whole family of solutions to the steady-state cooling flow 
equations derived in Nulsen (1986). It lies at one extreme 
with the following distribution at the other: 

/2 = ^^Po-'^e-(-/-)=-". (11) 



V2-a 



/2 gives a much more extended mass-deposition profile, M oc 
r"^ , at variance with the observations, but it is difficult to rule 
it out on purely theoretical grounds (it is possible that the 
comoving assumption may break down when gas phases of 
arbitrarily low density are included in the emulsion). 

In this paper we use fi and /2 to represent the extremes 
of the possible multiphase behaviour found in clusters. We 
emphasize once again that the cooling time in the bulk of 
a cluster is very long, so that one would not expect the 
steady-state cooling flow equations to give an accurate rep- 
resentation of the system. Nevertheless, the above distribu- 
tions provide a reasonable estimate of the expected range of 
behaviour. 



3 METHOD 

3.1 The cluster model 

We will suppose that the distributions of gas and mass in 
the cluster are determined via X-ray observations. The emis- 
sivity profile is often parameterized as 

/ 2\-3/3 

«(r)={o 1 + ^ , (12) 



where the core radius, a, and slope, /?, can be determined 
by fitting to the projected surface brightness profile. Jones 
& Forman (1984) give typical values a < 200/t~^kpc and 
/? ~ |. Direct deprojection of the surface brightness profiles 
(as described by Fabian et al. 1981) gives very similar re- 
sults. Thus the gas is much more extended than the galaxy 
distribution which declines as r~'^ in the outskirts of the 
cluster. 

Our first constraint on the system is therefore 

Ao < n'T'^ >= + ^) ' (13) 

where <> represents an average over all density phases, 
<A>= [ Afdp. (14) 



A second constraint is required in order to determine the 
density and temperature at all radii. In principle this could 
be provided by a direct measurement of the temperature 
profile, but unfortunately the observations are generally in- 
adequate to provide more than an overall spectrum for the 
whole cluster (spatially-resolved spectra are now becoming 
available from ASCA but even these are sufficient only to 
indicate the sense of the temperature variation, whether ap- 
proximately isothermal or declining with radius). We there- 
fore assume a polytropic relation. 



P x< p>'' 



(15) 



where 1 < 7 < 5/3. The lower limit corresponds to an 
isothermal cluster and the upper to one which is on the 
verge of convective instability. We ignore the effect of cool- 
ing flows which may give a temperature decline within the 
core of the cluster. 

The emission-weighted temperature averaged over the 
whole cluster is assumed to be determined by the X-ray 
observations. 



Tx 



J < n^T'^ > 47rr2dr 



(16) 



For the purposes of illustration we also define an emission 
weighted temperature as a function of radius. 



< T >e 



> 



(17) 



The distribution of gas in the cluster can be related to the 
total mass via the Equation of Hydrostatic Support which, 
for a spherically-symmetric system, has the following form. 



dP 
dr 



<P> 



GMtot 



(18) 



where G is the Gravitational constant and Mtot {t) is the 
total mass interior to radius r. 

Equations 13, 15, 16 and 18 together determine the dis- 
tribution of both the gas and the mass within the cluster. 
The single-phase equations are obtained by merely omitting 
the average over density phases. Before discussing the solu- 
tions, we will first convert the equations into dimensionless 
form. 



3.2 The dimensionless equations 

We define dimensionless variables (denoted by primes) as 



follows: 




r = 


1 

ar 


n = 


r]n' 


T = 


TxT' 


P = 


vkTxP' 


Mgas = 




Mtot = 





< n 1 > r dr 



With these definitions. Equation 16 becomes 

< n 1 > r dr . 

Choosing r] to satisfy = Ao»?^T^ gives 
< n 1 >= (1 -I- r ) ^. 



(19) 



(20) 



(21) 
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Figure 1. The temperature, density and cumulative gas and total 
mass profiles for a single-phase, isothermal cluster. 



Finally, we set //gas = 47ra iituhV l^tot = akTx I Gjimn ^ 
giving 



and 

ML 



< ra > r dr 



dP' 



(22) 



(23) 



< ra' > dr' ' 

The relative normalization of Mgas and Mtot is dependent 
upon fo, Ao and a whole host of observational parameters. 
Here we are concerned only with the difference between 
the normalization under the single-phase and multiphase as- 
sumptions. 

For clarity in the discussion that follows we will hence- 
force drop the primes in the above equations. Thus all vari- 
ables are to be considered dimensionless. 



4 RESULTS 

4.1 The isothermal case 

We start with the simpler case in which the icm is assumed 
to be isothermal (i.e. there is the same emission-weighted 
temperature at each radius). Figure 1 shows the variation of 
density, temperature, gas mass and total mass with radius 
for a single-phase icm. In this and the following figure we fix 
a = |- and /? = |: their precise values make little difference 
to our conclusions. We stress once again that the relative 
normalization of the two mass curves in our model is ar- 
bitrary. Thus, although the baryon fraction can be seen to 
increase as one moves out through the cluster, its magnitude 
cannot be determined from this plot. 

Now consider the multiphase case. As temperature is 
independent of radius, we have from Equation 20 



< n 1 



> = < n 1 > 



(24) 



(recall that we are omitting the primes for clarity of expres- 
sion). Hence 



1 — a ^ ^ 2 

< ra > = < ra 

or 

P =<n> Ct 
where 



> 



(25) 



(26) 



Table 1. Scaling factors for the gas mass, total mass and baryon 
fraction, relative to the single-phase case, for the multiphase vol- 
ume fraction distribution /i . 



a 






Cm 




(3-c.)l/ = 


2(2-a) 






2"/2(2-c.) 


(3-a) 


2(2 + <»)/2(2-a)2 





0.87 


1.33 


0.65 


1 

2 


0.89 


1.20 


0.74 



Table 2. Scaling factors for the gas mass, total mass and baryon 
fraction, relative to the single-phase case, for the multiphase vol- 
ume fraction distribution /2 . 



(2-a)2r(^ 



r 





2 



0.80 
0.82 



(2— )r(^) 

1.57 
1.37 



(2-a)f 



0.51 
0.60 



Ct 



< ra >< ra^ 



> 



(27) 



is a dimensionless constant which depends upon the func- 
tional form of the volume fraction distribution. 

From Equation 21, we see that the emissivity depends 
upon 



, 2 — a ^ 

< n 1 >= F < ra > = < ra > — ^ „ Cj. 



< ra >2 



Hence 

< ra >= C„(l + r^y 
where 



■3/3/2 



< ra 



< ra^ 



> C° 



< ra >< ra 



> 



a/2 



< ra2-« >(i + a)/2 



(28) 



(29) 



(30) 



The multiphase solution is therefore just a simple scaling 
of the single-phase one with the gas density being reduced 
by a factor C„. Similarly, from Equations 23 and 26, it is 
easy to see that the total mass is scaled upwards by a factor 
Ct relative to the single-phase case. The baryon fraction is 
therefore lowered by a factor 



Ca 



Ct ~ < ra2-« >(3+a)/2 



(31) 



Tables 1 and 2 give general expressions for the form-factors 
C„, Ct and Cm, plus numerical values for two different val- 
ues of a, for the two volume fraction distributions discussed 
in Section 2.2. The canonical form, /i, offers a reduction of 
26 percent to the inferred baryon fraction at high tempera- 
tures, rising to 35 percent in lower temperature clusters for 
which a lower value of a is more appropriate. The second 
form of the volume fraction, /2, that has a higher propor- 
tion of low-density phases, gives an even greater reduction 
of between 40 and 49 percent. 
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Figure 2. The temperature, density and cumulative gas and total 
mass profiles for a single-phase, polytropic cluster with 7 = 11/7. 



4.2 The polytropic case 

Next we investigate a more general temperature distribution 
of the form given by Equation 3. More specifically, we write 



P : 



K < n > 



1 



(32) 



where k is constant across the cluster. Increasing 7 results 
in a higher central temperature and hence a higher inferred 
gravitational potential. Thus the baryon fraction is lowered 
in the cluster core. At large radii, however, the reverse is 
true. Indeed, the condition 



7 



1 < 



2 



6/? ■ 



(33) 



must hold if the enclosed mass is not to decline (unphysi- 
cally) at large radii. Figure 2 shows the equivalent profiles 
to Figure 1 for the limiting case 7 = 11/7 Si 1.57 [a = |-, 
/? = I). The central baryon fraction is just one quarter of 
that in the isothermal case, but by r = 3 it is approximately 
equal and at larger radii it increases rapidly. Moving to a 
polytropic temperature law does not, therefore, help to al- 
leviate the baryon catastrophe. 

Although the formal derivation is harder than for the 
isothermal case, the form factors which correct the esti- 
mates of the gas and total mass are identical. This is be- 
cause the ratio of the 'dynamical' temperature, P / < n >, 
to the emission-weighted temperature depends only on the 
volume fraction distribution and not on the variation of P 
and < ra > across the cluster. 

To derive this result, consider first Equation 21. The 
left-hand-side can be written as 



< n 1 > = P < n > 



, 2 — a ^ 
a ^ " > ^ ^ 2-|-a(7-l) 

K < ra > ^ ' 

< n >2-« 



Hence 

<„>=A(l + r^)-^''/(^+'^(-'-^» 
where 



A 



< ra >" 



2-a \ 24-c(t-1) 



< ra2-« > K« 



(34) 



(35) 



(36) 



Contrast this with Equations 29 and 30 for the isothermal 
case. 



Equation 35 can now be substituted into the tempera- 
ture constraint. Equation 20, to yield after some manipula- 
tion 



7-1 



C'tL 



TJ-1 



where Ct is given by Equation 27 and 



I. 



Jr2(l -I- r2)-3/3(7+l + a(7-l))/(2-|-a(7-l)) 



(37) 



(38) 



is a constant independent of the volume fraction distribu- 
tion, /. At a given radius, the dependence of the total mass 
estimate on / comes from Equation 23, 

Aftot oc oc K < ra >''"V kA''~^ . (39) 

< ra > ' ' 

Hence from Equation 37 we see that the correction factor is 
equal to Ct, just as for the isothermal case. 

Finally, Equation 37 can be substituted back into Equa- 
tion 36 to eliminate k, giving 

<n>= C„I.r"l\l + ,2)-3/3/(2 + a(7-l))^ (40) 

where once again the multiphase correction factor, C„, is 
identical to the isothermal one given by Equation 30. 

In summary, moving to a polytropic temperature law 
makes little difference to our conclusions. The baryon frac- 
tion within the core of the cluster can be lowered, but at 
the expense of a higher fraction in the cluster as a whole. 
The correction factors which need to be applied to the mass 
estimates when moving from a single-phase to a multiphase 
icm are identical to those for the isothermal case. 



5 DISCUSSION 

We have shown that, if the icm is multiphase throughout, as 
we know it to be within the cooling flow at its core, then es- 
timates of the baryon fraction from X-ray observations need 
to be reduced by one quarter to one half. Taking a polytropic 
temperature distribution changes the estimate of the baryon 
fraction, making it lower in the cluster core (for 7 > 1) but 
substantially higher at large radii, however the correction 
factor when moving from a single-phase to a multiphase icm 
is identical to the isothermal case. 

In our model we have neglected a whole gamut of ob- 
servational details. In particular we considered the bolomet- 
ric flux emitted by the cooling gas. In reality this will be 
modified both by absorption, which tends to cut off the 
spectrum below a keV, and by the energy window of the 
detector. It is doubtful whether this would alter our conclu- 
sions greatly. The emission temperature in rich clusters is 
typically 3 x 10^-10** K which gives a flat spectrum from 
bremsstrahlung emission across the energy band of most 
imaging detectors. In a multiphase icm there is a mixture 
of gas at both higher and lower temperatures which will 
contribute a lesser or greater fraction of its emission in the 
observed energy band, respectively. As all density phases 
contribute in roughly equal measure to the observed emis- 
sion (the high-density phases are cooling most rapidly but 
occupy a much smaller volume) then the nett effect will be to 
leave the inferred bolometric emissivity almost unchanged. 

A correction factor of 26 or even 49 percent is not 
enough to resolve the baryon catastrophe all on its own. 
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Some fraction of the Universe could be made up of hot dark 
matter (hdm) that does not cluster efficiently into cluster 
halos, but the hdm mass-fraction cannot be greater than 
about 30 percent which is still not sufficient to save the idea 
of an Einstein-de Sitter Universe. 

Many simulations of cluster formation have been car- 
ried out with little evidence that baryons will concentrate 
preferentially in clusters. Both White et al. (f 993) and Babul 
& Katz (f993) looked specifically for this effect but found a 
baryon enhancement of only about 20 percent in the cluster 
as a whole. It should be pointed out, however, that there 
could be a lot of missing physics in these simulations, for 
example baryons could be congregated by inhomogeneous 
ionizing radiation fields in the early universe. 

Optical observations of cluster masses are prone to 
large errors, both through projection effects and because 
the three-dimensional shape of the velocity-dispersion ten- 
sor is unknown. However, X-ray mass determinations are 
much more reliable. Evrard (f994) suggests an accuracy of 
within ±50 percent, and Schindler (f995) an rms dispersion 
of f5 percent about the correct value, with no systematic 
under- or over-estimate. 

One further way of increasing cluster mass estimates, 
more in the spirit of this paper, is to modify the equation 
of hydrostatic support for the icm by adding extra terms to 
the pressure (e.g. Loeb & Mao f 994). The extra contribution 
from magnetic fields and turbulence in equipartition could 
triple the inferred masses, thus negating the Baryon Catas- 
trophe at a stroke. This hypothesis is supported by mod- 
elling of gravitational arcs which gives masses 2-2.5 times 
larger than those deduced from an unmodified equation of 
hydrostatic support (e.g. Miralda-Escude & Babul f995) al- 
though a correct modelling of shear may lower this estimate 
somewhat (Bartelmann, Steinmetz & Weiss f995). 

In a contrary vein, cold gas which has condensed from 
the icm will add to the observed baryon content of clusters 
and if it is coupled to the hot gas (a situation which has been 
termed 'mass loading') then it will lower the cluster mass 
estimate (by raising the mean density but not the pressure 
of the gas). We can only assume that this does not happen to 
a significant extent outside the cluster core, although there 
is evidence that it does do so within the cooling flow region 
(Daines, Fabian & Thomas f994). 

Our own personal view is that the intracluster medium 
is much more complex than most people have hitherto as- 
sumed and that there is sufficient uncertainty in its mod- 
elling so as to permit a critical density, Einstein-de Sitter 
Universe. 
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